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ABSTRACT 

In  this  thesis  several  incompressible  oscillatory  flow  and 
flutter  problems  were  investigated.  First,  a  previously 
developed  unsteady  panel  code  was  modified  so  that  systematic 
comparisons  with  Theodorsen's  classical  theory  could  be 
accomplished.  It  was  found  that  the  panel  code  is  in  excellent 
agreement  with  the  Theodorsen  results.  Second,  the  panel  code 
was  applied  to  the  analysis  of  bending- torsion  flutter.  Again, 
general  agreement  with  Theodorsen's  flutter  predictions  was 
obtained.  In  the  experimental  part  of  the  thesis  two  flow 
visualization  experiments  were  performed.  First,  the  vortical 
flow  patterns  generated  by  an  airfoil  executing  harmonic 
plunge  oscillations  were  visualized.  In  the  second  experiment, 
the  interference  effects  between  a  stationary  airfoil  and  a 
small  vane  executing  plunge  oscillations  were  explored. 
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I .  INTRODUCTION 

A.  GENERAL 

In  this  thesis,  several  numerical  methods  were  used  to 
analyze  the  flow  about  an  airfoil  performing  unsteady  motion 
in  an  inviscid  incompressible  fluid.  First,  the  unsteady 
motion  of  a  single  airfoil  was  analyzed  after  modifying  the 
U2DIIF  code  [ref.2].  The  primary  purpose  was  to  verify  the 
code  against  the  proven  theory  of  Theodorsen  for  analyzing  the 
phenomenon  of  flutter.  To  accomplish  this  the  U2DIIF  code  was 
modified  to  calculate  aerodynamic  values  over  a  range  of 
reduced  frequencies  and  then  apply  these  values  to  the 
flutter  analysis. 

Next,  the  propulsive  effects  of  a  plunging  airfoil  were 
verified  through  experimental  methods  using  a  low  speed 
plexiglas  wind  tunnel . 

Finally,  an  exploratory  test  was  conducted  in  the 
department's  smoke  tunnel  to  study  the  interaction  between  a 
plunging  airfoil  and  a  stationary  large  airfoil. 

B .  SCOPE 

Chapter  II  contains  the  modification  of  the  single  airfoil 
U2DIIF  code  into  the  code  UPOT.f  and  extensive  verification  of 
this  code  against  results  produced  by  Theodorsen.  Chapter  III 
describes  the  UPOT  code  and  explains  the  modifications  which 


were  added  to  solve  the  flutter  determinant.  In  chapter  IV  the 
flow  visualization  experiment  is  described  which  was  performed 
to  study  the  vortical  wake  patterns  produced  by  a  plunging 
airfoil .  In  chapter  V  a  second  experiment  is  described  which 
was  performed  to  explore  a  plunging  airfoil's  potential  for 
control  of  flow  separation. 


II.   SINGLE  AIRFOIL  ANALYSIS 

A.   U2DIFF  PANEL  CODE 

1 .  Geometry 

Figure  2.1  shows  a  representation  of  the  system  that 
is  analyzed  using  the  panel  code.  Shown  are  the  values  for  h 
(plunge)  and  a  (AOA) . 

2.  U2DIFF 

The  U2DIIF  code  was  developed  by  TENG  [ref  .2]  for  the 
study  of  unsteady  inviscid  and  incompressible  flow  over  a 
single  airfoil.  The  code  is  based  on  the  extension  of  the 
panel  method,  developed  by  Hess  &  Smith  [ref. 4]  for  steady 
potential  flow  problems,  to  include  the  unsteady  motion  of  the 
airfoil  that  is  continuously  shedding  vortices  into  the 
trailing  wake.  This  vortex  shedding  process  is  nonlinear  in 
that  the  wake  vortices  influence  the  flow  over  the  airfoil 
which  in  turn  alters  the  vortex  shedding  as  the  airfoil 
proceeds  in  time. 

The  non-linearity  of  the  unsteady  flow  makes  this 
problem  different  from  the  steady  flow  problem  which  requires 
only  simple  Gaussian  elimination.  Teng  developed  a  code  that 
used  an  iterative  type  of  solution.  Typical  program  output 
includes  the  airfoil  pressure  distribution,  force  and  moment 
coefficients,   and  the  trailing  vortex  wake  pattern.    No 
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Figure  2.1  Airfoil  Geometry 


attempt  is  made  here  to  reproduce  the  work  of  Teng  or  to 
explore  the  operation  of  the  U2DIIF  code,  but  the  reader  is 
encouraged  to  review  reference  2 . 

B.  PHASE  PROGRAM 

The  phase  program  was  put  through  some  verification  by 
Neace  [ref.9]  and  modified  slightly  in  order  to  present 
results  for  harmonic  motion.  The  code  PHV3 . f  (phaseshift)  was 
written  by  Neace  to  convert  the  time  dependent  output  of  lift 
and  moment  histories  to  harmonic  output  using  an  iterative 
curve  fit  algorithm: 

F{t)=Amp*Sin{<i>t+Q)  (2.1) 

where  Amp  =  amplitude  of  motion,  w  =  frequency,  and  0  =  phase 
angle  between  motion  and  the  aerodynamic  forces.  One  primary 
output  of  this  program  was  the  values  of  phaseshift  (0) 
between  the  AOA  and  coefficients  of  lift  (CL)  and  moment  (CM) 
for  the  pitching  airfoil  and  the  phaseshift  between  the 
plunge  value  (h/2b)  and  the  CL  and  CM  for  the  plunging  airfoil. 
The  other  output  was  the  amplitude  of  CL  and  CM  for  the 
pitching  or  plunging  case. 

C.  MODIFICATION  OF  U2DIIF  AND  PHASE  PROGRAM 

In  an  attempt  to  make  the  above  mentioned  codes  more  "user 
friendly",  the  two  codes  were  combined  into  a  single  code 
named  UPOT.f.   The  modification  involved  a  new  input  file 


called  UPOT.in  which  gives  the  user  several  options  of 
operation.  The  input  file  can  call  for  the  analysis  of  steady 
flow  only,  straight  and  modified  ramp  motion,  pitch 
oscillation,  plunge  oscillation,  and  the  capability  of 
performing  the  oscillation  analysis  over  a  series  of  reduced 
frequencies.  A  sample  input  file  is  shown  in  Figure  (2.2). 
1 .   Output 

Outputs  from  the  code  have  been  limited  to  reduce  the 
amount  of  computer  space  taken  up  by  the  code  operation.  A 
sample  calculation  was  run  using  the  input  from  Figure  (2.2) 
and  the  on-screen  output  is  shown  in  Figure  (2.3) .  Of  course, 
the  user  can  modify  the  output  portions  of  the  code  to 
minimize  output.  The  following  list  describes  the  input/output 
files  and  the  data  they  contain. 

a.  UPOT.IN:   The  input  file  figure  (2.2). 

b.  CL.d:  This  file  contains  the  various  AOA  values  along 
with  its  corresponding  CL  for  each  time  step. 

c.  CM.D:  This  file  contains  the  various  AOA  values  along 
with  its  corresponding  CM  for  each  time  step. 

d.  PHASE,  d:  This  file  contains  the  values  of  non 
dimensional  time  (t)  ,  AOA,  CL,  CM,  for  each  time  step. 

e.  FOR015.DAT:  This  file  contains  the  values  of  non 
dimensional  AOA,  curve  fit  for  CL,  curve  fit  for  CM, 

(used  in  the  phase  portion  of  program) . 

f .  CPSS.d:  This  file  contains  the  steady  state  pressure 
coefficient  for  the  mid  point  locations  of  all  the  air 
foil  panels. 

g.  CPU005.d:  This  file  contains  the  unsteady  pressure 
coefficient  for  the  mid  point  location  of  all  airfoil 


panels,  (in  this  case  the  values  are  for  an  AOA  equal 
to  5  degrees) . 

h.  PHZSWP.d:  This  file  contains  the  phase  information  of 
the  reduced  frequency  sweep  portion  for  the  program. 
The  file  contains  the  phase  angle  of  CL,  and  CM,  and 
the  amplitude  of  CL,  CM. 

i.  FLUTTER. IN:  This  file  contains  information  that  can 
be  used  to  solve  the  flutter  determinant.  It  contains 
K  ,  CL  Re,  CL  Im,  CM  Re,  CM  Im  for  the  pitch  or  plunge 
case. 


stdin 


Page  1 


AIRFOIL  TYPE 
NI.OWER  =  50  , 


:  NACA  0012  AIRFOIL 
NUPPER  =  50 


I  FLAG 

N LOWER 

NUPPER 

0 

50 

50 

AIRFOIL  TYPE 

IRAMP 

IOSCIL 

ALP  I      ALPMAX 

PIVOT 

0 

1 

-3.0       3.0 

0.37 

FREQ 

RFQSTP 

RFQFNL 

.68 

0.01 

.7 

IGUST 

UGUST 

VGUST 

0 

0. 

0. 

ITRANS 

DELHX 

DELHY  DELI   PHASE 

0 

.00 

.0     .0    0.00 

CYCLE 

vlTCYCLE 

TOL 

2 

60 

0.005 

naot  4 

naot  X 

aoa  values  multiplied 

by  1 C  ( integer  1 

2 

05  10 

20  25  39  50 

Comment  s  .  . . 

IRAMP   0:  -i/a  RFREQ  is  based  on  full  chord 

1 :  Straight  ramp 
2 :  Modi  f led  ramp 

IOSCIL  0:  r./a  RFREQ  is  based  on  full  chord 

1:  Sinusoidal  pitch,  motion  starts  at  run  Aoa 

ITRANS  0:  n/a 

1:  Trans lat lonal  harmonic  oscillation 

ALPI/ALPMAX    Minimum/MAX   AOA  in  degrees  for  IRAMP/ 1  TRANS/ IOSCIL 

MAX  does  not  apply  for  ITRANS 
PIVOT   Location  of  Elastic  Axis  as  a  fraction  of  full  chord 
FREQ    Initial  reduced  frequency  for  program 
RFQSTP  Reduced  freq  step  size  for  a  sweep  of  f req. ' s (enter  0.0  if  only  one  calcu 

lation  Is  desired.) 
RFQFNL  Final  freq  for  the  sweep 

DELHX   Translat ional  amount  In  the  chordwise  direction  (dist/full  chord) 
DELHY   Max  Translat lonal  amount  in  the  vertical  direction  (h/f ul lchord (b) ) 
DELI    Min  Translat ional  amount  in  the  vertical  direction  (h/b) 
CYCLE   :  I  of  cycles  for  oscillatory  motions 

-In  case  of  ramp,  cycle=l.S  denotes  airfoil  is  held 

at  max  aoa  for  the  duration  of  .5  cycle 
-For  steady  state  solution  set  it  to  0 

NTCYCLE:  I  of  time  steps  for  each  cycle 

CYCLE'NTCYCLE  is  limited  to  200  currently. 
TOL      Tolerance  for  convergence  of  the  unsteady  solut ion. (recommend  using  not 

less  than  .001) 
NAOT:  •  of  Input  aoa  for  cp  output 

-  angles  should  be  In  increasing  order, 

-  for  oscilatory  motions  angles  should  increase 
first,  then  decrease.  Decreasing  angles  are  for 
the  return  cycle.. 


Figure    2.2    UPOT.IN 


stdin 

Page  1 

4 

a:  rfo 

L  TYPE  :  NACA  0012  AIRFCIL 

NLOWER  =  50  ,  NUPPER  =  50 

IFLAG    (C:NACA,  1:INPUT) 

0 

NO. 

PANELS  UPPER  SURFACE      =    50 

NO. 

PANELS  LOWER  SURFACE      =    50 

7 
OSCILLATORY  MOTION,  IOSCIL 

INITIAL  ANGLE 

OF  ATTACK 

= 

■3.0000 

FINAL    ANGLE 

OF  ATTACK 

= 

3.0000 

REDUCED  FREQ. 

FOR  OSCIL 

= 

0. 6800 

REDUCED  FREQ. 

STEP 

0.0100 

FINAL  REDUCED 

FREQ. 

= 

0.7000 

PIVOT  POINT 

= 

0.3700 

TOTA 

L  1  OF  CYCLES 

= 

2.0C0C 

•  of 

TIME  STEPS  PER  CYCLE 

60 

TOLERANCE  FOR 

CONVERGENCE 

0.0050 

FREQ  SWEEP 

FREQ  = 

0.680000 

STEADY  FLOW  SOLUTION  AT  ALPHA  =   -: 

.000000 

1 

0.999507 

-1.116717 

-0.086820 

0.280753 

-0 

848C84 

2 

0.997535 

-1.108171 

-0.086820 

0.213074 

-0 

887088 

3 

0.993600 

-1.100064 

-0.086820 

0.168200 

-0 

912031 
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96 

2.5981 

14.4179 

1 

0.1625 

0.0005 

0.0322 

97 

2.4271 

14.5697 

0 

0.1461 

-0.0002 

0.0323 

98 

2.2294 

14.7215 

0 

0.1279 

-0.0008 

0.0320 

99 

2.0074 

14.8732 

0 

0.1081 

-0.0013 

0.0313 

100 

1.7634 

15.0250 

0 

0.0870 

-0.0016 

0.0303 

101 

1.5000 

15.1768 

0 

0.0647 

-0.0018 

0.0289 

102 

1.2202 

15.3285 

0 

0.0415 

-0.0018 

0.0272 

103 

0.9271 

15.4803 

0 

0.0178 

-0.0016 

0.0252 

104 

0.6237 

15.6321 

0 

-0.0063 

-0.0013 

0.0229 

105 

0.3136 

15.7838 

0 

-0.0305 

-0.0008 

0.0203 

106 

0.0000 

15.9356 

1 

-0.0546 

-0.0002 

0.0175 

107 

-0.3136 

16.0874 

0 

-0.0782 

0.0005 

0.0145 

108 

-0.6237 

16.2391 

0 

-0.1011 

0.0012 

0.0113 

109 

-0.9270 

16.3909 

0 

-0.1230 

0.0020 

0.0079 

110 

-1.2202 

16.5427 

0 

-0.1437 

0.0028 

0.0044 

111 

-1.5000 

16.6945 

1 

-0.1630 

0.0035 

0.0009 

112 

-1 .7634 

16.8462 

0 

-0.1807 

0.0042 

-0.0026 

113 

-2.0074 

16.9980 

0 

-0.1965 

0.0048 

-0.0061 

114 

-2.2294 

17.1498 

0 

-0.2103 

0.0052 

-0.0096 

115 

-2.4270 

17.3015 

0 

-0.2219 

0.0055 

-0.0130 

116 

-2.5981 

17.4533 

1 

-0.2312 

0.0056 

-0.0163 

117 

-2.7406 

17.6051 

0 

-0.2382 

0.0055 

-0.0194 

118 

-2.8532 

17.7568 

0 

-0.2427 

0.0053 

-0.0223 

119 

-2.9344 

17.9086 

0 

-0.2446 

0.0049 

-0.0250 

120 

-2.9836 

18.0604 

0 

-0.2441 

0.0044 

-0.0274 

PHASE  SHIFT  ANALYSIS 

FREQ  = 

0.6900000 

Figure    2.3g  UPOT   output 
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stdin 


Page  8 


AMPLITUDE;   clanp,  cmarrp 
PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


0.2300961  3.4413978E-02 

193.5791  -37.85255 

i.5437672E-03  9. 4169796E-02 

-0.2118947  -7.2855391E-C3 


Figure  2.3h  UPOT  output 
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2.   Difference  Between  UPOT  and  U2DIIF/Phase 

The  input  file  format  was  changed  along  with  the 
following: 


The  program  can  now  analyze  a  pitch,  plunge  or  ramp  motion 
that  starts  from  any  minimum  value  of  Alpha  or  plunge 
(h/2b)  .  Previously,  the  program  only  accepted  the  initial 
position  of  zero.  This  program  does  not  need  to  go 
through  the  origin. 


•  The  phase  portion  of  the  program  was  changed  to  curve  fit 

isine  function: 

F(  t)  =Amp*Cos  (co  t+4>)  (2.2) 


CL  and  CM  to  a  cosine  function: 


where,  Amp  =  amplitude  of  motion,  u  =  frequency  of  motion, 
<t>  =  phase  angle  between  motion  and  the  aerodynamic 
forces.  This  was  done  since  the  alpha  and  plunge  values 
were  allowed  to  start  from  a  new  zero  position. 

•  The  phase  portion  uses  the  middle  180  degrees  of  the  final 
360  degree  cycle  specified  in  the  UPOT. in  file.  This 
change  was  done  to  capture  an  all  positive  area  of  the 
cosine  curve  for  phasing  analysis.  The  program  integrates 
this  portion  of  the  cosine  curve,  and  for  proper  code 
operation  the  area  under  the  curve  must  be  kept  to  one 
sign.  If  the  areas  of  integration  were  chosen  to  include 
both  sides  of  the  axis,  then  the  code  would  produce  errors 
near  90  and  270  degrees. 

3.   UPOT  Verification 

The  code  UPOT  did  not  incorporate  any  drastic  changes 
to  the  prior  codes,  but  the  original  code  had  never  been 
extensively  compared  to  prior  theories  over  a  wide  range  of 
reduced  frequencies.  When  conducting  these  comparisons,  it  is 
easy  to  become  confused.  This  section  will  go  through  the 
comparisons  slowly  to  help  alleviate  that  problem. 
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a-       Kpanel     (*p)    vs  •    KTheodorson^Kt^  • 

The  equation  for  reduced  frequency  is: 


G)2£>     K  _  <jib 
U  c      U 


K=m?  K=^  (2.3) 


where:   w  =  frequency  of  oscillation  (rad/sec) 
b  =  half  chord  (units  to  match  U) 
U  =  free  stream  velocity  (units  to  match  b) . 
The  difference  between  X-    and  Kt   lies  in  the  fact  that  K— 
calls  for  the  full  chord  and  Kt  calls  for  the  half  chord, 
hence,  it  is  important  to  remember   that  K^   is  twice  Kt . 
b.      Aerodynamic  Forces 

The  aerodynamic  forces  problem  of  simple  harmonic 

motion  about  an  equilibrium  position  was  solved  theoretically 

by  Theodorsen  in  NACA  TR-496  [ref .10]  and  outlined  by  Fung  in 

[ref.5].   The  complex  equations  were  simplified  using  the 

simple  harmonic  motion  equation  and  resulted  in  the  following: 

L=7ipJb3a)2(Lh^+[La-(A+a)Lh]a+[Lp-(c-e)Lz]p)ei(wt+4,i)  (2.4) 


M=*p£4G)2  (  [Mh-  ( 1  +a)  Lh]  ^  + 
2  b 

[Ma-(A+a)  (La+Mh)  +(-|+a)2LA]a  (2.5) 

+  [Ma-(  1  +a)  LB-  (c-e)Mz+  (c-e)  (1+a)  Lz]  P)  ei(wt+v 


L,  M  are  the  lift  and  moment  per  unit  span  of  the 
airfoil  about  the  elastic  axis,  b,  h/b,  a  and  Oi    (radians)  ,  are 
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shown  in  Figure  (2.1).,  Lh,  ha,  Lp,  and  Ma  are  defined  by 
Scanlan  [ref . 6 , pp. 412 -424]  for  various  values  of  Kt  and  e. 
This  analysis  will  not  cover  airfoil  aileron  combinations. 
Therefore  (3   becomes  zero  and  equations  2.4  and  2.5  reduce  to: 

L=rcp£3w2(Lh^+[La-(A+a)Lja)ei<wt+^)       (2.6) 


M=npb*o>2  (  [Mh-  ( -|  +a)  Lh]  -|  + 
[MB-(-|+a)  (La+Mh)+<|+a)2Lja)ei(wt+^ 


(2.7) 


The  UPOT  panel  code  used  the  following  equations 
in  defining  lift  and  moment: 


C\  =  — ^-=^+jfei<wt+*L)    (fr^tan-1^       (2.8) 
L    2qb   v  L      L  WL  RL 


C«=-^=W+I^ei,WttW    <t>m=tan-^        (2.9) 
4gjb2  #„ 


where  RL  and  IL  are  the  real  and  imaginary  parts  of  CL,  and  RM 
and  IM  are  the  real  and  imaginary  parts  of  CM.  For  the  same 
conditions  the  lift  (L)  and  moment  (M)  should  be  the  same  for 
both  the  panel  code  and  Theodorsen.  This  fact  allows  for 
comparison  of  the  magnitude,  real,  imaginary  and  phase  of  lift 
and  moment . 

For  lift:  Lt  (eqn  2.6)  equals  L_  (eqn  2.8) 
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7iPjb3G)2  [Lhh/b+  (La-  (%+a)  Lh)  a]  ei(wt+^'  =2qbsJR27Ileu«'t^)      (2  .  10) 
After  canceling   ei(ut+*L): 


7iPjb3w2  [Lh(h/b)  +  (Le-  (V2+a)  Lh)  a]  =2qb^R2L  +  I2L         (2  .n) 

CL=^L2+JL2 


For  pitch   case,    h/b   =    0,    2.11   reduces   to: 

7iPJb3o)2  [  (La-  (V2+a)  Lh)  a]  =2g£Ca  (2  .12) 

Substitute   Kt    =   bw/U      for  a>2    and  q     =     &pU2      into 
equation  2.12   which  gives: 

2nqbK2  (La-  (Vfe+a)  Lh)  a=2gjbCLa  (2  .  13) 

After   cancelling  and  substituting  Kp  for  Kt: 


^[L„-(V2+a)LJa=CLa  «•«> 


This  relationship  can  be  further  broken  down  into  the  real  and 
imaginary  parts: 

nK2 
Imag:  £  [iLa-  (V2  +  a)  iLb]  <x=CLasin  (<J>L)    (2.15) 

Plunge  Case:   a=0  using  equation  2.11  gives: 
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nK2 
Real:  £  [La-(V2+a)Lja=CLacos(4)L)     (2.16) 


npb3u2Lh^-=2qbCLh  (2.17) 

The  panel  code  uses  h/2b  for  analysis  because  it  uses  full 
chord  vice  half  chord.  Therefore  equation  2.17  becomes: 

2npb3G>2Lh(-^-)=2qbCLb  (2.18) 


Substituting  as  before  for  or 


2nqbK2Lh2  ( -£- )  =2qbCLh  (2.19) 

£  J—' 


Cancel  and  substitute  Kp: 


r2 


(l^P)  (JL)l  =r  (2.20) 


This  can  also  be  broken  up  into  imaginary  and  real  parts  as 
before. 

MOMENT: 

Equating   equations   2.7   and   2.9    results    in: 

npb4u2  (  [Mh-  (V2+a)  Lh]  £+[Ma-  (1/2+a)  (La+Mh)  +  (V2+a)  2Lh]  a) 

b  an  (2.21) 


^qb^R'+I, 
For  pitch:      h/b   =    0 


21 


Mp=4gb2^JT^ei(wt+*'")  (2.22) 


C    =JR2+I2  (2.23) 

resulting    in: 

KPjb4o)2(  [Ma-(Vz+a)  (La+Mh)  +(V2+a)2Lh]a)=4qb2CMa  (2.24) 

After   substituting   and  cancelling: 


an  AT,2 


'-  [Ma-  (V2+a)  (Lu+Mh)  +  ( V2+a)  2LJ  =C„ 


'Mo 


REAL        Mh    =    fc 


8 


JMAG;      Mh    =    0 


/V  TT    Z^ 

5  [tfB-(%+a)  (L0+V2)+(V2+a)2Lj=CWacos(4)Ma)       <2-26) 


3  [iMa-  (%+a)  (iLa)+(V2+a)2iLj=CMasin(<t)Ma)      <2.27) 


8 

For  plunge,    a;   =    0 ,    equation   21    reduces    to 
REAL;      Mh    =    % 


-jKpi-j^)  [Mh-(V2+a)Lh]  =0^003  (t^)  (2.28) 
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IMAG:      Mh  =  0 

7  K2p{Yb]   (1/2+'a)Lh=C«hSin(<t>Mh)  (2.29) 

Comparisons  are  shown  for  various  cases  of  pitch  and  plunge. 
The  tables  include  pitch  values  of  1  (Tables  2.2-2.3)  and  6.7 
degrees  (Tables  2.4-2.5),  plunge  (h/2b)  values  of  .01  (Table 
2.8)  and  .0833  (Table  2.6-2.7).  The  graphs  include  1  degree 
pitch  (Figures  2.7-12,  2.13-20),  6.7  degree  pitch  (Figures 
2.21-28,2.29-36),  . 01  h/2b  plunge  (Figures  2 . 53- 56) ,  and. 0833 
h/2b  plunge  (Figures  2.37-44,  2.45-52). 

4.   Results 

The  tables  and  graphs  show  that  the  panel  code 
predicts  the  Theodorsen  results  accurately.  An  initial 
question  that  was  first  addressed  for  the  comparisons  was  how 
many  cycles  to  use  for  good  consistent  phase  results.  Initial 
runs  were  made  at  several  different  cycle  values  and  the 
results  are  shown  in  Table  2.1.  It  can  be  seen  that  the  panel 
code  predicts  Theodorsen' s  results,  using  a  cycle  number  of 
three.  Increasing  the  cycle  number  just  takes  more  computer 
time  and  only  results  in  marginal  increases  in  accuracy. 

The  most  glaring  difference  appears  in  the  IM  c^ 
comparisons  of  Figures  2.47-48.  The  panel  code  drops  off 
sharply  at  the  higher  end  of  the  reduced  frequency  spectrum. 
The  reason  for  this  is  believed  to  be  due  to  the  magnitude  of 
h/2b  chosen  for  the  comparison.   The  code  was  rerun  for  a 
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-  

Comparison  of  Phase  Calculations  Using  Various  Cycles. 

(Pitch,  6.7  deg .,  NACA  0007,  .37c,  50  panels  top  and  bottom) 

1 

Kp 

ft  Cycles 

CI  Phase  Angle 

Cm  Phase  Angle 

CI  Amplitude 

Cm  amplitude 

1.00 

1 

182.0537 

-54.409 

0.4884 

0.08937 

1.00 

2 

208.1592!                 -46.424 

0.5109224 

0.083932 

1.00 

3 

206.001 

-44.313 

0.51527 

0.083169 

%  Diff.  2/3 

1.05%  j                 4.76% 

0.84% 

0.92% 

1.00 

4 

204.9365 

-43.33887 

0.51668 

0.082951 

1.00 

5 

204.3955 

-42.868171      0.5172907 

0.082865 

1.00 

6 

204.0596 

-42.58497 

0.517598 

0.0828221 

1.00 

7 

203.8291 

-42.3916 

0.517776 

0.082792 

1.00 

8 

203.6031 

-42.25684 

0.51789 

0.082784 

%  Diff.  7/8 

0.11% 

0.32%  |            0.02% 

0.01%: 

3.60 

2 

264.44052 

59.79002 

1.0931 

0.212011 

3.60 

3 

261.9737 

-58.29977 

1.09954 

0.21936 

%  Diff.  2/3 

094% 

2.56% 

0.59% 

3.35% 

3.60 

4 

260.17877 

57.4404 

1.10271 

0.21903 

3.60 

5 

259.2686 

56.79395|        1.104505 

0.218845 

3.60 

6 

258.5186 

56.3252 

1.105599 

0.218739 

3.60 

7 

257.9483 

56.0098        1.106306 

0.218673 

3.60 

8 

257.6163 

•55.380567 

1.106764 

0.218632 

%  Diff.  7/8 

0.13% 

1.14%  |            0.04%  i       0.02% 

TABLE  2.1  PHASE  CALCULATION  VS  CYCLE  NUMBER 

series  of  h/2b  values  and  the  percent  difference  for  the  panel 
code  to  Theodorsen  was  plotted  in  Figure  2.4.  This  chart  shows 
that  the  h/2b  value  chosen  has  a  tremendous  impact  on  the  code 
results.  An  h/2b  value  of  .01  gave  an  acceptable  error  of  10% 
at  IL  =  8.  Runs  were  completed  with  a  value  of  .01  h/2b  and 
the  favorable  results  are  shown  in  Figures  2.53-56. 
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Real  Part  of  CI  for  Panel  and  Theodorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure  2.5  1  Degree  pitch  CT  Re 
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Real  Part  of  CI  for  Panel  and  Theodorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc  per  cycle) 
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Figure  2.6  1  degree  pitch  CL  Re 
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Imag.  Part  of  CI  for  Panoi  and  Theordorsen  (pitch,  1.0  deg,  .37c, 

N ACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure  2.7  1  Degree  pitch  CL  Im 
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Imag.  Part  of  CI  for  Panel  and  Theordorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Mag.  of  CI  for  Panel  and  Theordorsen  (pitch,  1.0  dag,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc  per  cycle) 
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Mag.  of  CI  for  Panel  and  Theordorsen  (pitch,  1.0  dag,  .37c, 
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NACA0007, 
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and  bottom,  3c 

:ycles  of  65  ca 
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Phase  of  CI  for  Panel  and  Theordorsen  (pitch,  1.0  deg,  .37c, 
NACA0007, 75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure  2.11  1  Degree  pitch  CL  Phase 
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Phase  of  CI  for  Pane!  and  Theordorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Cm  Real  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure  2.13  1  Degree  pitch  CM  Re 
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Cm  Real  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1.0  deg,  .37c, 
NACA0007,  75  panels  top  and  bottom,  3cycles  of  65  calc.  per  cycle) 
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Cm  Imag  vs  Kp  for  Panel  and  Theordorsen  (pitch,  1 .0 

deg,  .37c,  NACA0007,  75  panels  top  and  bottom, 

3cycles  of  65  calc.  per  cycle) 
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Cm  Imag  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1 .0 

deg,  .37c,  NACA0007,  75  panels  top  and  bottom, 

3cycles  of  65  calc.  per  cycle) 
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Cm  Magnitude  vs  Kp  for  Panel  and  Theordorsen 

(pitch,  1.0  deg,  .37c,  NACA0007,  75  panels  top  and 

bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure   2.17    1   Degree  pitch  CM  Magnitude 
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Cm  Magnitude  vs  Kp  for  Panel  and  Theodorsen  (pitch, 

1 .0  deg,  .37c,  NACA0007,  75  panels  top  and 

bottom,  3cycles  of  65  calc.  per  cycle) 
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Figure   2.18    1   Degree  pitch  CM  Magnitude 
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Cm  Phase  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1 .0 

deg,  .37c,  NACA0007,  75  panels  top  and  bottom, 

3cycles  of  65  calc.  per  cycle) 
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Figure   2.19    1   Degree  pitch  CM  Phase 
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Cm  Phase  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1 .0 

deg,  .37c,  NACA0007,  75  panels  top  and  bottom, 

3cycles  of  65  calc.  per  cycle) 
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Figure  2.20  1  Degree  pitch  CM  Phase 
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CI  Real  vs  Kp  Panel  Code  (Pitch.  6.7  deg,  .37c,  NACA  0007,  50  panels 
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CI  Real  vs  Kp  Panel  Code  (Pitch,  6.7  deg,  .37,  NACA  0007, 75  panels 
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CI  Imag.  vs  Kp  Panel  Code  (pitch,  6.7  deg,  .37c,  NACA  0007, 75  panel 

top  and  bottom) 
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CI  Imag.  vs  Kp  Panel  Code  (pitch.  6.7  deg,  .37c,  NACA  0007,  75  panel 
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CI  Magnitude  vs  Kp  Panel  Code  (pitch,  6.7  deg,  .37c,  NACA  0007.  50 
panels  top  and  bottom) 
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CI  Phase  vs  Kp  (pitch,  6.7  deg,  .37c,  NACA  0007, 75  panel  top  and 

bottom) 
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Figure   2.27    6.7    degrees  pitch  CL  phase 
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Figure  2.31  6.7  Degrees  pitch  CM  Im 
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Figure  2.32  6.7  Degrees  pitch  CM  Im 
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Figure  2.33  6.7  Degrees  pitch  CM  magnitude 
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Figure   2.34    6.7   Degrees  pitch  CM  magnitude 
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Cm  Phase  vs  Kp  Panel  Code  (pitch,  6.7deg,  .37c, 
NACA  0007,  50  panels  top  and  bottom) 
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Figure   2.35    6.7   Degrees  pitch  CM  phase 
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Figure    2.3  8    Plunge   h/2b=.0833    CL  Re 
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Figure  2.41  Plunge  h/2b=.0833  CL  Magnitude 
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Figure   2.42    Plunge  h/2b=.0833      CL  magnitude 
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Phase  of  CI  for  Panel  and  Theodorsen  (plunge,  .0833 

h/2b,  .37c,  NACA0007,75  panels  top  and  bottom, 

3cycles  of  65  calculations) 
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Figure   2.43    Plunge   h/2b=.0833   CL  phase 
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Figure   2.44    Plunge  h/2b=.0833    CL  phase 
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TABLE    2.7    PLUNGE   h/2b=.0833    CM   COMPARISON 
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Real  Part  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.45  Plunge  h/2b=.0833  CM  Re 
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Real  Part  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007.75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.46  Plunge  h/2b=.0833  CM  Re 
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Imag.  Part  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.47  Plunge  h/2b=.0833  CM  Im 
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Imag.  Part  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007.75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.48  Plunge  h/2b=.0833  CM  Im 
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Magnitude  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure   2.49    Plunge  h/2b=.0833    CM     Magnitude 
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Magnitude  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.50  Plunge  h/2b=.0833  CM  Magnitude 
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Phase  of  CM  for  Panel  and  Theodorsen  (plunge, 

.0833  h/2b,  .37c,  NACA0007.75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure   2.51   Plunge   h/2b=.0833    CM     Phase 
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Phase  of  CM  for  Panel  and  Theodorsen  (plunge. 

.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure  2.52  Plunge  h/2b=.0833  CM  Phase 
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TABLE  2.8  PLUNGE  h/2b= . 01  CM   COMPARISON 
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Real  Part  of  CM  for  Panel  and  Theodorsen  (plunge, 

.01    h/2b,  .37c,  NACA0007,100  panels  top  and 

bottom,  3cycles  of  65  calculations) 
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Figure   2.53    Plunge   h/2b=.01   CM     Re 
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Imag  Part  of  CM  for  Panel  and  Theodorsen  (plunge,  .01 

h/2b,  .37c,  NACA0007,100  panels  top  and  bottom, 

3cycles  of  65  calculations) 
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Figure   2.54    Plunge  h/2b=.01   CM     Im 
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Magnitude  of  CM  for  Panel  and  Theodorsen  (plunge,  .01 

h/2b,  .37c,  NACA0007,100  panels  top  and  bottom. 

3cycles  of  65  calculations) 
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Figure  2.55  Plunge  h/2b=.01  CM  Magnitude 
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Phase  of  CM  for  Panel  and  Theodorsen  (plunge,  .01 

h/2b,  .37c,  NACA0007,100  panels  top  and  bottom, 

3cycles  of  65  calculations) 
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Figure  2.56  Plunge  h/2b=.01  CM  Phase 
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Comparison  of  CM  Imag  vs  Kp  for  Various  Values  of  Pitch 
(NACA0007,100  panels,  1,3and  6.7  deg.  pitch,3cyc65)(0c.) 
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Figure  2.57  Pitch  Oc.  1.0, 3.0, and  6.7  degrees  vs  Theodorsen 
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'Comparison  of  Imaginary  Dimensionless  Aerodynamic  Coefficient  (Im 
My=8*lm  Cmalpha|(alpha*pi*Kp'2)(NACA0007,100  panels,  1,3.6.7  deg. 

pitch,3cyc65,0c.) 
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Figure  2.58  Dimensionless  Aerodynamic  Coefficient  for 
1.0f3.0fand  6.7  degrees 


86 


III.   FLUTTER  DETERMINANT 

The  proven  accuracy  of  the  UPOT  Code  enabled  it  to  be  used 
for  the  solution  of  the   flutter  determinant. 

A.   FLUTTER  THEORY 

In  order  to  analyze  the  phenomenon  of  flutter,  it  is 
necessary  to  obtain  the  equations  of  motion  of  the  system.  To 
simplify  the  problem  the  assumption  is  made  that  the  actual 
motion  of  the  system  can  be  considered  a  combination  of 
fundamental  wing  bending,  and  fundamental  wing  torsion.  The 
system  can  then  be  replaced  by  an  equivalent  system  containing 
an  airfoil  section  of  unit  span  restrained  by  springs  against 
independent  vertical  motion  (bending) , and  torsion  as 
illustrated  in  Figure  3.1.  This  paper  will  not  consider  the 
aileron  hinge  case  so  £  and  c  are  set  equal  to  zero.  According 
to  the  class  notes  of  M.  Platzer  [ref.l]  the  formulation 
proceeds  as  follows: 

Consider  the  balance  of  the  elastic,  inertial  and 
aerodynamic  forces  on  a  mass  element: 

•  Total  Inertial  force:   -/  dm(h"+ra" ) =- (Mh"+Saa" ) 
Mass :   M=/  dm 

Static  Moment  about  the  elastic  axis:  Sa  =  /  rdm 

•  The  moments  about  the  elastic  axis  are: 
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-!r(h"+ra.")dm  =  -  {Iaa"+Sah")  (3.1) 

Mass  moment  of  inertia  about  elastic  axis 
Ia  =  jr2dm 


Elastic  restoring  forces  are:   -hCh     ~aCa 
The  Equations  of  motion  therefore  become: 


h"M+a"Sa+hCh  =  L 


a."la+h"Sa+a.Ca  =  M 


(3.2) 


Where:  Ca  =  Torsional  stiffness  of  the  wing 

Ch  =  Stiffness  of  the  wing  in  translation  (plunge) 
M  =  Mass  of  the  wing  per  unit  span 
These  equations  can  be  written  in  a  different  way  by 
expressing  the  spring  constants  in  terms  of  the  natural 
frequencies.  Consider  the  airfoil  to  be  so  restrained  that 
only  one  degree  of  freedom  is  permitted.  The  equations  of 
motion  become: 


C 
Mh"+hCh   =  0  so  that      o>„  =  x  — 

h       \    M 


(3.3) 


Iacrl0+aC0  =  0      so  that      u^  = 


Hence:         Ch   =   Mcoh2         Ca   =    Iawa2 


N  J- 


The  small  structural  damping  of  metal  aircraft  may  be 
approximated  by  a  force  that  opposes  the  motion  and  is  in 
phase  with  the  velocity.  One  assumes  therefore  that  the 

88 


magnitude  of  the  damping  is  proportional  to  the  elastic 
restoring  force.  Since  the  motion  of  the  airfoil  is  harmonic 
at  the  critical  flutter  condition,  the  structural  damping  can 
be  accounted  for  by  replacing  the  terms: 

hCh  with  hCh(l  +  igh) 

oiCa  with  aCa(l  +  iga) 
Where  gh  and  ga  are  damping  constants  multiplied  by  i  to  ensure 
that   the  damping  force  is   in  phase  with  the  velocities   in  the 
simple   harmonic   motion. 
From  equation   3.3    with: 

h(t)    =   heiot      and     a(t)    =  aeiot 
We   have : 

h"    =   o2heiut   and     a"    =    -w2aeiot 
And  the   equations   of  motion  become: 

eiot(-fa2hM-<a2aSa+hCh)    =  L  (3.4) 


eiut(-w2aJB-a)2h5tt  +  aC(,)    =M  (3.5) 


The  equations  for  the  aerodynamic  forces  were  given  by 
Fung  [ref.5]  and  are  shown  here: 

L=7ip£3w2(L^+[La-(-^+a)Lja)eiwt         (3.6) 

Equating  equation  3.4  to  3.6  and  3.5  to  3.7  yields: 
Substituting  into  equation  3.8  and  3.9  for  Ch  and  Ca  and 
using  the  following  dimensional  terms: 
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M=npb*u2(  [Mb-(±+a)Lb]  |  + 
[Mtf-(l+a)  (La^)+(A+a)2Lja)eiwt 


(3.7) 


(-w2M2-a)2a5/hq)    =rcp£3w2  (Lh^  +  [La- (  — +a)  Lh]  a)      (3.8) 


(-o2aIs-u2ii5B+«cJ  =np2?4a>2(  [Mh-(±+a)Lh]  ^  + 
[Ma-  ( |  +a)  (La+MA)  +  ( 1  +a)  2LJ  a ) 


(3.9) 


H  = 


M 


rcpir 


5j 

"a  Mb 


x  =  a 


N 


j. 


Mb: 


(3.10) 


The  equations  simplify  and  after  bringing  all  terms  to  the 
left  result  in: 


A^+Ba=0 

D 

D^+Ea=0 

D 


(3.11) 


This  is  a  homogeneous  equation  whose  solution  is  obtained 
if  the  flutter  determinant  is  zero. 


A  B 

D  E 


=   0 


(3.12) 


Where 
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A  =  VL  [i  -  (  _^  )  2  (  _^Z  )  2  ( i  +  igr  )  ]  +Lh 
(0  0)a  fl  fl 

B=liX+L  -Lh(1/2  +  a) 

H*     «        «        ft  (3.13) 

D=jixa+1/2-Lh(1/2+a) 

E=|ira2[l-(  —  )2(l+iga)]  -1/2(V2+a)  +Ma-La(V2+a)  +Lh(V2+a)2 

CO 

/x  is  the  ratio  of  the  mass  of  the  wing  to  the  mass  of  a 
cylinder  of  air  of  a  diameter  equal  to  the  chord  of  the  wing. 


u> 


and  ii>h  are  the  natural  angular  frequency  (rad/sec)  of 
torsional  vibration  around  "a"  (elastic  axis)  and  the  natural 
frequency  in  deflection,  respectively.  xa  is  the  location  of 
center  of  gravity  of  the  wing  measured  from  a.  to  is  the 
circular  frequency  of  wing  vibration. 

The  relationships  between  the  code  and  Theodorsen  derived 
earlier  in  Chapter  II  can  be  used  here  to  simplify  the 
equations:   (note:   no  damping  in  this  case  ga=gh=0) 

For  A:  manipulating  equation  2.20, 


r      .  ZCLh 


h     -„2<h_,  (3.14) 

2b 


nK2p(^L) 


resulting  in 


A=lUi-(.^M-gV]+  ZC».  (3  15) 


*^p(^r) 


2b 
For  B:  manipulating  equation  2.14, 
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4  a 


La-Lh{V*+a)=-^f-  (3.16) 

nKp(x 

resulting   in: 

B=HXa  + f-  (3.17) 

nKpd 

For  D:  manipulating  equation  2.28 


resulting  in 


D=UX°* TfhT  (3-19) 


For  E:  manipulating  equation  2.24 


8  C, 


-1/2(1/2+a)  +Ma-La(V2+a)+L/:(V2+a)2  = ^-       (3.20) 

an*/ 

resulting  in: 

i?=^ra2[l-(-^)2]+-^-  (3.21) 

W        CC7liCp2 

The  determinant  is  expanded  to  AE-BD=0,  and  the  real  and 
imaginary  parts  are  set  equal  to  zero.  Substituting  (coa/oj)  = 
X  and  solving  the  real  (2  roots)  and  imaginary  (1  root) 
equations  for  values  of  X  corresponding  to  each  reduced 
frequency  value.   These  X  values  can  be  plotted  as  SQRT(X) 


92 


against  K  and  any  intersections  of  real  and  imaginary  parts 
signify  a  flutter  point. 
Knowing  that  : 


jr=i£».      and  ^=(  — )        <3-22> 

p       U  a) 


solve    for  Uerfticll 


2jbo>„ 

^critie^  — ~=  (3.23) 


Pv 


which  is  the  critical  flutter  speed. 

B.   UPOTFLUT  CODE 

1.   FORMULATION  AND  INPUT 

The  equations  derived  in  the  flutter  theory  section 
above  were  programmed  into  a  FORTRAN  subroutine  and  attached 
to  the  UPOT.f  code.  The  UPOT  code  was  modified  first  to 
enable  it  to  conduct  a  frequency  sweep  of  pitch  and  plunge 
simultaneously.  The  resulting  frequency  sweep  pitch  and 
plunge  array  data  is  then  sent  to  the  flutter  subroutine  which 
provides  the  values  of  SQRT(x)  and  Kp  for  plotting.  The 
program  also  gives  a  best  guess  for  the  Ucritical  based  on  the 
difference  betweeen  the  real  and  imaginary  SQRT  (X)  values. 
The  input  file  UPOTFLUT.  IN  is  very  similar  to  the  regular 
UPOT. IN  file  with  the  addition  of  actual  physical  properties 
of  the  system  being  analyzed.    The  user  should  start  the 
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analysis  in  the  pitch  mode  first  (IOSCIL  =  1,  ITRANS  =0)  to 
ensure  complete  coverage  of  all  frequencies  of  interest.  The 
following  relations  were  taken  from  NACA  TR-685  [ref.8]  and 
should  prove  helpful  in  determining  the  physical  properties 
needed  for  program  operation. 

k    =   mass  ratio  =  irpb2/M 

k    =  l//z 

xa  =  Sa/Mb 

ra2    =  Ia/Mb2 

2 .   OUTPUT 

Outputs  from  the  code  have  been  limited  to  reduce  the 
amount  of  computer  space  taken  up  by  the  code  operation.  A 
sample  input  and  output  file  are  contained  in  Figure  3.2. 
The  following  list  describes  the  input/output  files  and  the 
data  they  contain: 


a.  UPOTFLUT.IN  The  input  file  figure  3.2a 

b .  CL .  d  Same  as  UPOT . f  output 
Same  as  UPOT.f  output 
Same  as  UPOT.f  output 
Same  as  UPOT.f  output 
Same  as  UPOT.f  output 
This  file  contains  Kp,  <£L,  <£M,  CLa,  CMa 

h.  PLHZSWP.d   This  file  contains  K^ ,  0L0M ,  C^ ,  C^ 

i.  PITCH,  in    This  file  contains  YL^    ,  CLa  Re,  CLa  ImlClto  Re, 
CMaN  Im 


CM.d 
PHASE. d 
CPSS.d 
CPU005.d 
PHZSWP.d 
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'imik 


m 


Wing  Bending 
Spring 

Wing  lorsion  Spring 
at  Rotation  Point 


Midchord 


■i^E^ 


Aileron  Torsion 
Spring  at  Aileron  Hinge 


b  =  semichord  (ft.) 

cb  =  distance  between  midchord  and  aileron  hinge,  positive  if  aft  of 
midchord  (ft.) 

eb  =  distance  between  midchord  and  aileron  leading  edge,  positive 
-  aft  of  midchord  (ft.) 

ab  =  distance  between  rotation  point  (elastic  axis)  and  midchord, 
positive  if  aft  of  midchord  (ft.) 

h  =  bending  deflection  of  rotation  point  (elastic  axis),  positive  down- 
ward (ft.) 

t  a  =  angular  deflection  about  rotation  point  (elastic  axis),  positive 
for  leading  edge  up  (radians) 

/9  =  angular  deflection  of  aileron  about  aileron  hinge  relative  to  wing 
chord,  positive  for  aileron  leading  edge  up  (radians) 


Figure  3.1  Simplified  System  Geometry 
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stdln 


Page  t 


AIRFOIL  TYPE  :  NACA  0012  AIRFOIL 
NLOWER  =  50  ,  NUPPER  «  50 


ALPMAX 
3.0 


PIVOT 
0.3 


IFI.AG   NLOWER   NUPPER 
0         50    50 
AIRFOIL  TYPE 

7 
IRAMP   IOSCIL    ALPI 

0       1      -3.0 
FREO    RFQSTP   RFOFNL 

.85      .01      .95 
IGUST   UGUST   VGUST 

0      0.       0. 
ITRANS   DELIIX   DELHY   DELI   PHASE 
0     0.00    .0833   -.0833   0.00 
CYCLE   NTCYCLF      TOL 

3      65       0.005 
naot  (   naot  X  aoa  values  multiplied  by  10  (Integer) 

2       05  10   20  25  39  50 
Semi-chord   Walpha    Wh      Mass 

6  90       72.0      .53789 

Ialpha      Salpha   Density 
4.84102     .645468    .002378 
Comments . . . 


IRAMP   0:  n/a 


RFREO  is  based  on  lul 1  chord 


1:  Straight  ramp 
2:  Modified  ramp 


IOSCIL  0:  n/a  RFREQ  is  based  on  full  chord 

1:  Sinusoidal  pitch,  motion  starts  at  min  Aoa 

ITRANS  0:  n/a 

1:  Trans lat lonal  harmonic  oscillation 

CYCLE   :  I  of  cycles  for  oscillatory  motions 

-In  case  of  ramp,  cycle=1.5  denotes  airfoil  Is  held 

at  max  aoa  for  the  duration  of  .5  cycle 
-For  steady  state  solution  set  It  to  0 

NTCYCLE:  •  of  time  steps  for  each  cycle 

CYCLE'NTCYCLE  is  limited  to  200  currently. 

NAOT:  I  of  input  aoa  for  cp  output 

-  angles  should  be  in  Increasing  order, 

-  for  oscllatory  motions  angles  should  Increase 
first,  then  decrease.  Decreasing  angles  are  for 
the  return  cycle.. 

SEMI-C5H0RD  Half  Chord  in  feet. 

Walpha, Wh,  uncoupled  natural  frequencies  of  the  system  In  question. 
Walpha  is  pitch  and  Wh  Is  plunge(HZ). 

Mass        specific  mass  of  the  system  in  slugs/foot  of  span 

Ialpha      Moment  of  Inertia  of  system  about  the  elastic  axls(a) 

per  unit  span  length. 
Salpha      Static  moment  of  wing-aileron  per  unit  span  length 

Density     Mass  of  air  per  unit  of  volumefslugs  per  ft"3) 


Figure   3.2a  UPOTFLUT.In   example   input   file 
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stdin 


AIRFOIL  TYPE  :  NACA  0012  AIRFOIL 
NLOWER  =  50  ,  NUPPER  >  50 


OSCILLATORY  MOTION,  IOSCIL 


FREQ  SWEEP 

FREO  *    0.850000 


w   0.8500000 
kp  =   0.8500000 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.8500000 


1  f  req 


AMPLITUDE;   clamp,  cmatnp  : 
loscll  =  Oltrans 

PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


0.2316691 


2.9221054E-07 


202.4033       -63.70797 
1.5930884F.-03   0.1051438 
-O.H30808      -1.1134184E-02 


w   0.8500000 
kp=   0.8500000 


PHASE  SHIFT  ANALYSIS 
FREO  »    0.8500000 


1  f  req 


AMPLITUDE;   clamp,  cmamp  : 
loscll  =  Oltrans 

PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREO  SWEEP 

FREO  =    0.860000 


0.2779089 


7.2117507E-02 


1 


270.5030       37.10253 

-6.2051453E-03  -0.4095396 

0.1050769      -0.1181067 


v      0.8600000 
kp=   0.8600000 


PHASE  SHIFT  ANALYSIS 
FREQ  =    0.8600000 

lfreq  2 


AMPLITUDE;   clamp,  cmamp  :   0.2317267 
.ioscll  =  Oltrans  = 


2.9490557E-02 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


202.8760      -63.85251 
1.5921656E-03   0.1050829 
-0.1398697      -1.1383206E-02 


w   0.8600000 
kp  =   0.8600000 


PHASE  SHIFT  ANALYSIS 
FREQ  =    0.8600000 


lfreq 


AMPLITUDE;   clamp,  cmamp  : 
ioscll  =  Oltrans 

PHASE;      clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREQ  SWEEP 

FREQ  =    0.870000 


0.2804866 


2.2494521E-02 


1 


270.7843       36.79393 

-6.3188463E-03  -0.4170439 

0.1047920     -0.1205979 


w   0.8700000 


PHASE  SHIFT  ANALYSIS 
FREQ  '         0.8700000 


Figure  3.2b  UPOTFLUT  example  output  file 
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stdiri 


kp  =      0.8700000 


I  f  req 


AMPLITUDE;   clamp,  cmamp  : 
loscll  =  Oltrans 

PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


0.2317936 


2.9765518E-02 


203.3604       -64.02047 
1.5893428E-03   0.1048966 
-0.1365460      -1.1639614E-02 


w   0.8700000 
kp  =   0.8700000 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.8700000 

If req  3 


AMPLITUDE;   clamp,  cmamp  :   0.2830637 
loscll  =  Oltrans  =  1 


2.2874046E-02 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREO  SWEEP 

FREO  =    0.880000 


271.0655       36.47753 

-6.4334869E-03  -0.4246101 

0.1045149      -0.1231114 


w   0.8800000 
kp-   0.8800000 


PHASE  SHIFT  ANALYSIS 
FREO  »    0.8800000 

If req  4 


AMPLITUDE;   clamp,  cmamp  :   0.2318744       3 .00371 80E-02 
loscll  =  Oltrans  =  0 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


203.8330       -64.15525 
1.5880425E-O3   0.1048108 
-0.1335113      -1.1894441E-02 


w   0.8800000 
kp  =   0.8800000 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.8800000 


1  f  req 


AMPLITUDE;   clamp,  cmamp  : 
loscll  =  Oltrans 

PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREO  SWEEP 

FREO  ■    0.890000 


0.2856959 


2.3258235E-02 


1 


271.3546        36.17284 

-6.5488103E-03  -0.4322215 

0.1042214      -0.1256711 


w   0.8900000 
kp=   0.8900000 


PHASE  SHIFT  ANALYSIS 
FREQ  =    0.8900000 

If req  5 


AMPLITUDE; 
losctl  = 


clamp,  cmamp  :   0.2319494 
Oltrans  = 


3.0311935E-02 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


204.3115       -64.29196 
1.5859993E-03   0.1046759 
-0.1304959      -1.2153637E-02 


w   0.8900000 
kp=   0.8900000 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.8900000 


AMPLITUDE; 


If req  5 

clamp,  cmamp  :   0.2883391 


2.3644408E-O2 


Figure   3.2c   UPOTFLUT   example   output   file 
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stain 


ioscll  »  Oltrans 

PHASE;       clp,    crop 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREQ  SWEEP 

FREO  =    0.900000 


1 
271.6436        35.87401 
-6.6649993E-03  -0.4398900 
0.1039310      -0.1282581 


w   0.9000000 
kp=   0.9000000 


PHASE  SHIFT  ANALYSIS 
FREO  =■    0.9000000 

1 f rcq  6 


3.0592278E-02 


AMPLITUDE;   clamp,  cmamp  :  0.2320280 

Ioscll  =            Oltrans  =            0 

PHASE;       clp,    crop    :  204.7803       -64.43259 

AVERAGE  DRAG,  TOTAL  DRAG  :  1 . 582881 1E-03   0.1044/01 

ETAS,  WBAR                 :  -0.1274617      -1  .  24 1 84  82F.-07 


w   0.9000000 
kp=*   0.9000000 


PHASE  SHIFT  ANALYSIS 
FRF.Q  '  0.9000000 

lfreq  6 


AMPLITUDE; 
loscil  - 


clamp,  cmamp  :   0.2909827 

Oltrans  =  1 


2.4039967E-02 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREQ  SWEEP 

FREO  =    0.910000 


271.9288        35.58691 

-6.7821071E-03  -0.4476191 

0.1036481      -0.1308680 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.9)00000 


w   0.9100000 
kp  =   0.9100000 


1  f  req 


AMPLITUDE;   clamp,  cmamp  : 
loscil  =  Oltrans 

PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


0.2321466 


3.0865142E-02 


205.2529       -64.61227 
1.580O473E-O3   0.1042831 
-0.1245367      -1.2687406E-02 


w   0.9100000 
kp=   0.9100000 


FHASE  SHIFT  ANALYSIS 
FREO  *         0.9100000 

lfreq  7 


AMPLITUDE; 
loscil  « 


clamp,  cmamp  :   0.2936268 

Oltrans  =  1 


2.4440434E-02 


PHASE;       clp,    cmp 
AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


FREO  SWEEP 

FREO  »    0.920000 


272.2159        35.29980 

-6.9O00O78E-O3  -0.4554005 

0.1033707      -0.1335003 


w   0.9200000 
kp=   0.9200000 


PHASE  SHIFT  ANALYSIS 
FREO  =    0.9200000 

lfreq  8 


Page  3 


Figure   3. 2d  UPOTPLUT   example  output   file 
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AMPLIT'IDF."   cl-imp,  cmemp  : 

load!  =  Oltians 

PHASE;  dp,  cinp  : 
AVERAGE  DRAG,  TOTAI  DRAG  : 
FTA5,  WRAR  : 


i.7i77'5i  i.iiia697F-o? 

o 

205.1373  -64.49424 
1.5782372F-01     0.1011637 

-0.1718161  -1  .29191HF-0? 


PHASE    SMIPT    ANALYSIS 
TPFO    -  0.9700000 


v   O.97OOO0O 
kp-   0.9700000 


1  I  rr^ 


ArlPLITI'DF;   clomp,  cmamp 
loscll  =  Oltron 

PHASE;      clp,    cmp 
AVERAGE  DRAG,  TOTAI  DRAG 
FTA5,  WPAR 


FREO  swrrp 

rpro  -   0.910000 


0.7967705 


7.1R1158'F-02 


2 '2.5011       35.01268 

-7.0186830F-O1  -0. 4432331 

0.1010»B'      -0.1161516 


».   0.9799999 
Irp-   0.9711199 


PHASE  SHIFT  ANALYSIS 
FRFO  -    0.  9799999 


1  f  r*»l 


1.11171IOF-0? 


AMPLIT'IDF;       clomp,    cramp    :       0.7171518 

loscll     «  Olttonj    -  0 

PHASE;  rip,         rmp  204.2080  -41.82713 

AVERAGE  DRAG,    TOTAL   DRAG    :         1 . 5 755624F-03      0.10398)1 

ETAS.    HBAR  :       -0.1191*57  -1  .3771407F.-07 


PHASE  SIUFT  ANALYSIS 
IRPO  •    0.  9799911 


»   0. 9799999 
kp-   0.9799999 


AMP  LIT!  IDF 
losdl    -  Oltra 

PHASE;  clp,  •  mp 

AVERAGE    DRAG,     TOTAI.    DRAG 
F.TAS,    WRAR 


FREQ    SWEFP 

rpp.o  ■       0.910000 


I  fret]  9 

clump,    cmemp    I      0.7  9891 M 


7.5750117F-02 


2'?. 7901  31.72948 

-7.1382108F.-O3    -0.4711219 
0.1078376  -0.1388117 


-   0.9199999 
kp-   0.9199919 


PHASE  SHtTT  ANALYSIS 
PRFO  -    0. 9199199 


I  fieri 


in 


AMFLITUDF:       c  1  .-imp,    crimp    :  0.7374174                  1 .  I  6190  I6P.-0? 

loscll     »                               Oltrans  -                               0 

PHASE;               dp,         rvp  204.4814               -41.94821 

AVERAGE    DRAG,     TOTAI.    DRAG    :  1  . 5  715122r-0  J       0.1037211 

FTAS,     WRAR                                          :  -0.1164158               -1 . 14969H9E-07 


w      0.9'99999 
kp-       01194999 


PHASE    SHII I    ANALrSl 
FRFO    "  0.9199999 


1  f  f  rfj 


AMFLlTUDFr   clomp,  cmamp  : 
loicl]  •  01 1  fans 


in 
0.  1015610 


7.56S9967F-0? 


PHASE;       clp,    emp 
AVERAGE  DRAG.  T01AI  DRAG 
FTAS,  WRAR 

Number  of  Kp  Values  ■ 

Plunge  Value/Full  ChcrH  (li/7h) 

Alpha 

PIVOT  POINT  (•>  Of  Plastic  Aula   . 

Hall  Chord  (b) 

Alpha   uncoupled  Nat.    Iieq.  • 

Plunge   Uncouple^    Nar .    Freq 

I    alpha 

S   alpha  * 

rlaa* 

Alt    Penalty 

Man    Rat  lo 

Diet.     To    Ming    CG    Aft    of    Elal.    Aula 


273.0792  J4. 14823 

7.258U37r-03    -0.4190)51 
0.1075746  -0.1415101 

10 
0.0833 
3.0000 
-0.1000 
(.0000 
90.0000 
77.0000 

4.8410 
0.6455 
0.5379 
0.0074 

2.0000 
0.2OO0 


Kp  crlt  -  0. 9000000 

D1FT  ■  1.150)2)3F-nl 

SORTX  =  0.95)9428 

U  crlt  175; .484 


Figure   3.2e  UPOTFLUT   example   output    file 
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j.  PLUNGE. in   This  file  contains  K^,    CLh  RE,  C^  IM,  C^  RE, 

k.  FLUTPLOT.d  This  file  contains  Kp,SQRT(x)  Re,  SQRT(X)  RE, 
SQRT(X)  IM 

3 .   VALIDATION 

The  program  was  tested  against  some  sample  cases  to 
check  for  code  validity.  The  first  case  was  taken  from 
reference  6  example  #1,  p.  236.  Figures  3.3  and  3.4  show 
plots  of  the  FLUTPOT.d  file.  Figure  3.3  shows  the  initial 
look  over  a  wide  range  of  K_  and  after  finding  the  approximate 
flutter  location  Figure  3.4  shows  a  closer  look  at  the  K^ 
range  of  interest.  This  example  calculated  a  Ucritical  of 
161.985  ft/sec.  which  compares  favorably  to  the  example  value 
of  162  ft/sec.  The  next  example  was  taken  from  NACA  TR-685 
[ref.8]  case  #1  p.  8.  Figures  3.5  and  3.6  again  show  the 
initial  and  final  looks  for  this  analysis.  The  example  called 
for  a  Ucritical  of  567  miles/hr  and  the  program  returned  a 
value  of  570  miles/hr.  Next,  the  code  was  tested  over  a  range 
of  coh/coa  ratios  as  done  in  NACA  TR-685,  pll.,  graph  I-A(a)  . 
Figure  3.7  shows  the  comparison  between  the  two  methods. 
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Figure  3.3  Initial  look  at  flutter,  example  1 
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Figure  3.4  Final  look  at  flutter,  example  1 
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Figure  3.5  Initial  look  at  flutter,  example  2 
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Figure  3.6  Final  look  at  flutter,  example  2 
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Mechanism  of  Flutter   The  effects  ot  changing  the  ratio  of  plunge  and  pitch  natural  frequencies. 

Panel  code  vs  Theo.  (NACA0007,  BOpaneto.  3deg.  h/2b  =  .0833,  3cyc66calc,  ) 


Wh/Walpha 


0.0111 


0.2 


0.4 


0.6 


0.6 


0.7 


0  8 


(Mae»=  .63788  skips/ft.  Ialpha»  4.84102,  Sjlpna-  646468.b*=  ft,  std  day  air  density) 


Panel  Ucrit. (ft/sec) 


896.6013 


925  6568 


960.3143 


994.4232 


1036.167 


1107.237 


1262  889 


Panel  1-Ucrit/bWalphal 


Theo  (Ucrit/bWalpha) 


1.66 


-1.713994 


_-1.J7 

•1.72 


1.7783598 


-1.84162 


1.9188 


2.060438 


-1.8 


■I.J 


-2.04 


%Ditt  wrl  Theo. 


2.36% 


0.36% 


1.20% 


308% 


6  94% 


-2.22 


7.64% 


2.32016! 


-2.42 


4.13% 


Mechanism  of  Flutter 
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Figure   3.7    uh/ua   Calculations 
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IV.   FLOW  VISUALIZATION  EXPERIMENT 

A.   INTRODUCTION 

The  purpose  of  this  experiment  was  to  document  the 
production  of  thrust  by  a  plunging  airfoil.  This  was  a 
preliminary  experiment  to  better  understand  the  vortex  pattern 
produced  by  a  plunging  airfoil,  and  to  examine  the  production 
of  thrust  using  smoke  flow  visualization  techniques. 

An  explanation  of  what  constitutes  a  propulsive  vortical 
signature  along  with  smoke  flow  visualization  of  the 
propulsive  vortical  patterns  is  given  in  Reference  7.  In  this 
reference,  the  explanation  is  given  by  contrasting  the 
vortical  pattern  produced  by  a  cylinder  (drag)  with  the 
vortical  pattern  produced  by  a  plunging  airfoil  (thrust) .  The 
cylinder  produced  a  vortical  sheet  where  the  top  row  of 
vortices  rotated  clockwise  and  the  bottom  row  of  vortices 
rotated  counterclockwise.  This  pattern  induces  a  velocity 
component  in  the  upstream  direction  (Biot-Savart  law) .  In 
contrast,  the  plunging  airfoil  produced  a  clockwise  rotating 
vortex  sheet  on  the  bottom  row.  This  pattern  induces  a 
velocity  component  in  the  downstream  direction.  Reproduction 
of  the  flow  visualization  data  from  Reference  7  is  shown  in 
Figure  4.1. 
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Figure  4.1  View  of  flow  over  cylinder  (top)  and  plunging 
airfoil  (bottom)  [Ref.  7] 
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B .  THEORY 

A  comparison  was  done  using  the  incompressible  panel  code, 
U2DIIF.  The  purpose  of  this  study  was  to  examine  the  vortical 
pattern  produced  by  the  panel  code,  and  determine  if  the 
vortical  signature  matched  experimental  results.  The  input  to 
the  panel  code  was  set  up  to  best  match  the  conditions  of  the 
experiment  described  in  the  next  section.  The  panel  code  was 
run  using  a  plunge  amplitude,  h/2b  equal  to  .1977,  a  reduced 
frequency  of  1.8  and  a  zero  mean  AOA.  The  results  of  the 
vortical  pattern  are  shown  in  Figure  4.2.  Aside  from  the 
starting  vortex,  this  is  clearly  a  thrust  producing  vortical 
sheet.  Furthermore,  the  vortical  pattern  is  similar  to  that 
produced  by  the  experiment  shown  in  Figures  4.10  and  4.11. 

C.  EXPERIMENTAL  SETUP 
1.   Plunging  Airfoil 

The  plunging  airfoil  used  in  this  paper  was  originally 
a  wing  taken  from  the  rotor  of  a  model  helicopter.  The  wing 
was  attached  to  a  MB250  Shaker  Table  as  shown  in  Figure  4.3 
The  wing  was  made  from  a  NACA0007  airfoil  section  and 
consisted  of  a  2.45"  chord  and  a  22"  span.  The  wing  was  built 
from  a  foam  core  and  finished  with  a  layer  of  graphite  epoxy 
composite  for  added  fatigue  strength.  The  airfoil's  drive 
mechanism  was  a  MB  250  Shaker  Table  capable  of  1"  total 
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Figure  4.2  Wake  pattern  produced  by  U2DIIF  code 
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Figure  4.3  Shaker  table  with  wing 
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deflection.  The  Shaker  Table  was  limited  only  by  resonance 
frequencies  of  the  wing  which  occurred  around  20  Hz  or  1200 
rpm. 

2.   WIND  TUNNEL 

The  wind  tunnel  used  in  this  experiment  was  a  very  low 
speed,  low  turbulence  smoke  tunnel.  It  is  made  of  plexiglass 
walls  and  a  contraction  ratio  of  2.8:1.  The  motor  provides 
wind  tunnel  velocities  between  0  and  10  feet  per  second  (fps)  . 
The  smoke  was  created  using  a  Rosco  smoke  generator  and  piped 
into  the  tunnel  in  the  test  section  using  a  small  seven  tube 
smoke  rake  constructed  for  this  experiment.  Figure  4.4  is  a 
photograph  of  the  wind  tunnel  and  smoke  rake  used  in  this 
experiment . 

D.   TEST  PROCEDURE 

Testing  was  conducted  in  the  low  speed  smoke  tunnel  under 
several  different  conditions.  The  speeds  of  the  tunnel  were 
approximately  1.04  fps,  1.47  fps,  and  1.56  fps  (measured 
visually) .  These  low  speeds  allowed  good  pictures  and  the 
ability  to  get  higher  reduced  frequencies  without  calling  for 
too  high  a  load  on  the  wing.  The  actual  plunging  harmonic 
frequencies  ranged  from  1  to  15  Hz  and  amplitudes  from  1/16" 
to  1"  peak  to  peak.  The  tunnel  was  initially  turned  off  and 
the  Shaker  Table  turned  on  with  stagnant  smoke  in  the  tunnel . 
The  purpose  was  to  see  if  the 
plunging  airfoil  would  draw  the  smoke  through  the  tunnel  like 
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Figure  4.4  Tunnel  and  smoke  rake 
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a  fan,   thus  showing  the  production  of  thrust  by  the  plunging 
airfoil . 

Photos  were  taken  using  a  Nikon  35mm  camera  and  Kodak 
TMAX-400  ASA  black  and  white  film.  The  shutter  speed  was  set 
to  1/125  seconds  with  an  aperture  setting  of  4 . 0  for  the  light 
conditions.  Film  developing  time  was  optimized  at  9  minutes  at 
75  degrees  F. 

E.   RESULTS  AND  DISCUSSION 

The  result  for  the  tunnel  off  condition  flow  visualization 
experiment  was  as  expected.  The  wing  in  fact  accelerated  the 
smoke  in  its  vicinity. 

The  result  of  the  additional  rake  flow  visualization 
experiments  are  shown  in  Figures  4.5-4.14.  Figure  4.5  shows 
the  stationary  airfoil  at  zero  degree  AOA.  The  Reynolds  number 
(based  on  airfoil  chord)  is  10,000.  It  can  be  seen  that  the 
airfoil  produces  a  small  wake  with  the  boundary  layer  mostly 
attached.  Figures  4.6  through  4.14  show  the  vortical  wake 
flow  patterns  produced  by  plunge  oscillations  at  various 
frequencies  as  indicated.  Most  of  these  pictures  reveal  the 
propulsive  vortical  street  pattern  discovered  in  Reference  7. 
Previous  experiments  by  Neace,  [ref .9]  found  that  the  tunnel 
was  too  small  for  the  airfoil  size  used,  but  the  airfoil  size 
for  the  present  experiment  seemed  to  be  optimum,  as  seen  by 
the  long  trail  of  vortices.  The  vortical  patterns  show  that 
the  bottom  vortex  is  rotating  clockwise,  and  the  top  vortex  is 
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Figure  4.5  Steady  airfoil  1.56  ft/s 


Figure    4.6    Kp    =    1.8008,    h/2b=,1977,     1.56    ft/i 
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Figure   4.7    Kp    =    2.467,    h/2b=. 10204,    1.56    ft/: 
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"igure    4.8    Kp    =    2.467,    h/2b=.1913,     1.56    ft/s 
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Figure  4.9  Kp  =  4.112,  h/2b=. 14031,  1.56  ft/s 
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'igure  4.10  K^    =  6.167,  h/2b=. 05102,  1.56  ft/s 
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'igure    4.11    Kp    =    6.753,    h/2b=.1084,    1.56    ft/s 
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Figure  4.12  Kp  =  7.401,  h/2b= . 03826,  1.56  ft/s 
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Figure  4.13  Kp  =8.223#  h/2b= . 01913 ,  1.56  ft/s 
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Figure    4.14,    K^=12.335,    h/2b=. 01775,     1.56    ft/s 
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rotating  counterclockwise,  which  is  a  thrust  producing 
vortical  sheet.  It  can  be  seen  in  the  pictures  that  the 
frequency  greatly  affects  the  vortical  strength  (size) . 
Increasing  the  frequency  leads  to  an  increase  in  wake 
vorticity . 
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V.  LIFT  ENHANCEMENT  PRODUCED  BY  A  PLUNGING  AIRFOIL 

A.  THEORY 

Chapter  IV  demonstrated  the  propulsive  capability  of  a 
plunging  airfoil.  The  production  of  thrust  implies  the 
generation  of  a  jet  flow  which,  in  turn,  may  be  utilized  as  a 
boundary  layer  control  device.  Therefore,  an  additional  test 
was  conducted  in  the  NPS  smoke  tunnel  in  order  to  explore  the 
feasibility  of  this  concept. 

B .  SETUP 

The  same  NACA0007  plunging  wing  was  used  as  in  chapter  IV 
with  a  different  driving  mechanism.  The  wing  was  mounted  to 
an  ELECTRO- SEIS  Model  113  Shaker  Table  by  APS  Dynamics,  Inc. 
The  shaker  was  located  below  the  test  section  of  the  NPS  smoke 
tunnel  (Figure  5.1).  The  plunging  airfoil  was  mounted  to 
struts  at  both  ends  to  prevent  excessive  bending  while 
plunging.  The  large  airfoil  is  a  cambered  profile  taken  from 
the  rotor  of  a  full  size  helicopter  (13"  chord,  2"  Thickness, 
modified  NACA  airfoil,  Reynolds  Number  of  52,000) .  The  wing  is 
suspended  from  the  tunnel  ceiling  as  shown  in  the  flow 
pictures.  The  design  allowed  for  full  movement  of  the  big  wing 
to  position  it  in  the  vicinity  of  the  plunging  airfoil. 
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C.   WIND  TUNNEL 

This  study  used  the  Naval  Postgraduate  School's  flow 
visualization  wind  tunnel.  The  tunnel  is  an  open- circuit  one, 
with  air  entering  an  inlet  that  measures  4.5  m  X  4.5  m 
(15'X15').  As  the  air  enters  the  tunnel,  it  passes  through  a 
7. 5 -cm  long  honeycomb.  A  9:1  ratio  square  contraction  cone 
directs  the  flow  into  a  test  section  that  is  1.5  m  X  1.5  m 
(5'X5'),  and  6.7  m  (22')  long,  as  seen  in  Figure  5.2.  The 
flow  is  then  exhausted  into  the  atmosphere  through  a  fan, 
which  uses  variable  pitch  blades  to  control  the  speed  of  the 
flow.  The  speed  control  toggle  switch  is  located  right  below 
the  red  and  green  on/off  switch  located  in  the  left  side  of 
the  tunnel  control  room.  The  tunnel  speed  was  determined  using 
a  digital  manometer  which  was  verified  for  accuracy  (Figure 
5.3)  . 

An  observation  booth  is  located  on  the  side  of  the  tunnel . 
A  glass  window,  1.6  m  X  1.1  m  (5.2'  X  3.4'),  provides  the 
primary  viewing  area  from  the  observation  room  and  a  second 
one,  0.4  m  X  1.23  m  (1.33'X4')#  is  located  in  the  tunnel's 
roof.  The  main  viewing  window  had  sufficient  area  for  most  of 
the  photography,  with  the  top  window  used  for  illumination.  A 
circular  turntable  was  located  on  the  floor  of  the  test 
section  [ref.ll]  which  allowed  for  easy  access  to  the  shaker 
table.  The  walls  and  floor  of  the  test  section  were  flat  black 
for  low  light  reflectivity. 
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Figure  5.1  Shaker  table  setup  below  tunnel 
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Figure  5.2  Tunnel  Layout 
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Figure  5.3  Tunnel  speed  check 
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D.  SMOKE  GENERATION 

The  smoke  was  generated  in  the  Rosco  smoke/fog  machine. 
Many  different  smoke  injection  techniques  were  tried  but  with 
less  than  satisfactory  results.  Smoke  rakes  were  first  tried 
outside  the  tunnel  with  the  tube  number  varying  from  2  to  3  0 
tubes.  The  tubes  were  inserted  in  the  honeycomb  and  also 
separated  different  distances  from  the  inlet  of  the  tunnel. 
The  rake  was  also  tried  inside  the  test  section  with  very  bad 
results  (smoke  dispersed  immediately) .  Problems  ranged  from 
lack  of  smoke  and  turbulence  when  enough  smoke  was  present. 
The  Rosco  machine  at  its  lowest  setting  was  producing  a  very 
high  smoke  volume  and  whenever  the  flow  was  restricted  by  a 
smoke  rake  the  smoke  production  went  way  down.  The  final 
technique  adopted  was  very  simple.  The  smoke  output  was  sent 
directly  from  the  machine  to  a  1"  nozzle  which  was  manually 
waved  at  the  entrance  of  the  tunnel  to  make  a  steady  cloud. 
The  steady  cloud  was  gradually  pulled  into  the  tunnel, 
producing  a  thick  smoke  sheet  in  the  test  section. 

E .  PHOTOGRAPHY 

Photos  were  taken  using  a  Nikon  35mm  camera  and  Kodak 
TMAX-400  ASA  black  and  white  film.  The  film  speed  was  set  to 
1/250  seconds  with  an  aperture  setting  of  4.0  for  the  light 
conditions.  Film  developing  time  was  optimized  at  9  minutes  at 
75  degrees  F. 
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F.  EXPERIMENTAL  PROCEDURES 

The  first  step  involved  a  look  at  the  large  airfoil  to 
verify  normal  flow  patterns  (Figure  5.4  and  5.5)  and  find  the 
AOA  for  initial  trailing  edge  separation  (Figure  5.6).  Next, 
the  plunging  airfoil  was  placed  in  the  tunnel  by  itself  and  a 
run  was  made  to  verify  the  propulsive  capability  in  the  larger 
tunnel  at  higher  speeds.  As  seen  in  Figure  5.7,  the  airfoil 
produced  a  drag  vortical  flow.  Figure  5.8  and  5.9  shows  the 
propulsive  pattern  of  the  propulsive  airfoil.  Finally,  the  two 
airfoils  were  placed  in  close  relative  position  to  see  the 
interference  effect  between  the  two  airfoils. 

G.  RESULTS  AND  DISCUSSION 

Several  airfoil  position  combinations  were  studied,  as 
shown  in  Figures  5.10  through  5.15.  Figure  5.10  and  5.11  show 
the  plunging  airfoil  located  at  approximately  .65  chord  of  the 
large  airfoil.  Figure  5.12  through  5.15  show  the  plunging 
airfoil  located  at  approximately  .75  chord. 

The  differences  between  the  plunging  on  and  off  condition 
were  not  easy  to  see  with  the  eye  but  pictures  indeed  showed 
some  differences  between  the  two  conditions.  A  shortcoming  of 
this  experiment  was  the  inability  of  the  plunging  airfoil  to 
run  parallel  with  the  large  airfoil.  Additionally,  sizing  of 
and  relative  positioning  of  the  two  airfoils  was  not  optimized 
to  give  best  results.  The  two  airfoils  were  chosen  from  the 
resources  available  and  time  constraints  prevented  a  more 


132 


Figure  5.4  Large  airfoil  at  zero  AOA 
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Figure    5 


7TT^gT"airfoil    at    10    degrees   AOA 
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Figure  5.6  Large  airfoil  at  12  degrees  AOA  (stall) 
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Figure  5.7  Plunge  airfoil  steady 
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Figure  5.8  Plunging  airfoil  K=l*71 
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Figure  5.9  Plunging  airfoil  1^=3.42 
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Figure  5.10  Large  airfoil  at  12  degrees  AOA,  steady  plunging 
airfoil  position  1 
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Figure  5.11  Large  airfoil  at  12  degrees  AOA# 
airfoil   K=3.42,  position  1 


plunging 


.40 


Figure  5.12  Large  airfoil  at  12  degrees  AOA,   plunging 
airfoil  steady,  position  2 
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Figure  5.13  Large  airfoil  at  12  degrees  AOA, 
airfoil   K=3.42,  position  2 
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plunging 


Figure  5.14  Large  airfoil  at  14  degrees  AOA,   plunging 
airfoil  steady,  position  2 
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Figure  5.15  Large  airfoil  at  14  degrees  AOA#   plunging 
airfoil  K=3.42,  position  2 
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detailed  investigation  of  the  interference  effects  between  the 
two  airfoils. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  SINGLE  AIRFOIL  ANALYSIS 

The  modified  version  of  U2DIIF  (UPOT)  can  perform 
aerodynamic  calculations  over  any  range  of  reduced 
frequencies.  The  nonlinear  theory  presented  here  for  harmonic 
motion,  and  the  phase  relationships  that  exist  between  the 
airfoil  motion  and  the  aerodynamic  forces  have  been 
extensively  verified  by  comparison  with  Theodorsen's  linear 
theory.  Furthermore,  this  panel  code  was  applied  to  the 
analysis  of  incompressible  bending- torsion  airfoil  flutter. 
Again,  excellent  agreement  with  the  classical  Theodorsen 
analysis  was  obtained. 

Access  to  faster  computational  means  is  recommended  to 
shorten  the  time  needed  to  predict  the  flutter  points.  The 
code  should  be  modified  to  incorporate  three-dimensional 
calculations  which  would  help  solve  more  difficult  flutter 
problems. 

B.  FLOW  VISUALIZATION  EXPERIMENTS 

The  flow  visualization  experiment  successfully  showed  the 
development  of  thrust  produced  by  a  plunging  airfoil.  The 
enhanced  lift  experiment,  on  the  other  hand,  was  not  a 
complete  success.  The  smoke  visualization  presented 
difficulties  that  were  not  satisfactorily  overcome.  As  a 
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result,  the  pictures  taken  were  somewhat  inconclusive. 
Futhermore ,  the  angle  of  attack  of  the  oscillating  airfoil 
could  not  be  changed  thus  making  it  difficult  to  achieve  a 
flow  condition  conducive  to  lift  enhancement. 

It  is  recommended  that  further  experiments  be  conducted  in 
the  low  speed  smoke  tunnel  with  a  shaker  table  capable  of 
moving  an  airfoil  at  harmonic  frequencies  near  40  HZ. 
Additionally,  the  airfoil  must  be  modified  to  allow  change  of 
AOA.  Finally,  the  Rosco  smoke  machine  output  volume  must  be 
modified  to  permit  much  lower  smoke  output.  This  final  point 
proved  to  be  the  single  largest  detriment  to  the  visualization 
experiment . 
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